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(a)  Objectivei: 

The  baaic  objective  of  this  research  effort  has  been  to  develop  both  theory  and  algo¬ 
rithms  capable  of  providing  realistic  control  systems  for  physical  plants  of  interest  to  the  Air 
Force  which  are  jnodeled  as  distributed  parameter  systems,  i.e.,  infinite  dimensional  linear 
systems.  Applications  include  large  space  structures  and  precision  pointing  of  optical  and 
radar  systems. 

(b)  Accomplishments  and  Progress:  Executive  Summary 

In  this  section  ,  we  provide  a  summary  of  the  progress  and  accomplishments  under 
this  project.  More  detailed  discussions  are  included  as  Appendix  1. 

In  developing  feedback  control  systems  for  infinite  dimensional  linear  systems,  the 
major  stumbling  block  has  been  to  perform  the  inner-outer  factorization  (lOP):  one  must 
compute  an  all-pass  (inner)  factor  and  a  minimum  phase  (outer)  factor  in  order  to  solve  the 
control  problem. 

We  have  analyzed  the  lOF  and  have  developed  numerical  algorithms  to  accomplish  the 
factorization  and  have  studied  the  relationship  to  the  “well  posedness”  of  the  control  problem. 
We  have  also  been  able  to  relate  the  well-posedness  to  the  mixed  sensitivity  problem:  weU- 
posedness  implies  that  the  solution  to  the  mixed  sensitivity  problem  is  also  well  posed. 
Consequently,  we  now  have  a  clear  connection  between  the  lOF  and  solving  the  H"‘  mixed 
sensitivity  problem.  Conversely,  these  results  help  to  characterize  when  a  problem  is 
ill-posed,  i.e.,  when,  for  example,  an  all-pass  system  is  approximated  by  a  minimum  phase 
factor. 

The  criterion  we  have  developed  is  a  norm  bound  between  the  system  and  ones 
knowledge  of  the  system.  Let  {/„}  C  be  the  sequence  of  approximations  and  /  C  Haa 
the  true  system.  Then,  if 

7-00  '  1  +  f’  ' 

the  problem  is  well  posed,  and  the  sequence  of  outer  factorizations  converges  to  the  outer 
factor  of  the  true  system.  In  addition,  if  we  also  have  uniform  convergence,  /n  — *  f,  we 
also  get  the  inner  factorization  of  /n  to  converge  to  f. 

The  relationship  developed  between  lOF  and  well  posedness  also  contributes  to  the 
mixed  sensitivity  problem:  well-posedness  implies  the  solution  to  the  mixed  sensitivity  prob¬ 
lem  is  also  well  posed.  Additional  results  relate  to  developing  stabilizing  compensators  and 
the  interpolation  of  compensators.  (See  Appendix  1.) 
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We  have  also  extended  the  solution  of  single-input  single-output  infinite  dimensional 
systems  to  holdTor  a-wider  class  of  plants,  and  have  also  made  extensions  for  a  specific  class 
of  multi-input  multi-output  problems.  In  addition,  we  have  been  able  to  obtain  the  lOF  for 
a  specific  class  of  irrational  matrices.  FWther  details  are  given  in  Appendix  1. 
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Appendix  lA 


Research  Summary:  Katherine  M.  Crow 

This  research  examines  the  practical  implementation  of  direct  solution  of  /f*  optimal 
control  problems  for  distributed  parameter  systems.  We  have  outlined  a  procedure  for  nu¬ 
merical  computation  of  the  optimal  performance  and  optimal  compensator  starting  from  a 
general  plant  which  may  be  irrational  and  unstable.  In  fact,  this  procedure  can  be  applied 
to  a  plant  described  only  by  numerical  frequency  response  data.  In  addition,  we  have  speci¬ 
fied  conditions  under  which  the  problem  is  well-posed  in  terms  of  the  continuity  of  the 
performance  criterion  under  small  perturbations  of  the  problem  data.  These  results  relate 
to  the  numerical  conditioning  of  the  computational  procedure,  and  therefore  are  of  direct 
practical  significance.  In  order  to  motivate  this  work,  we  have  been  using  several  distributed 
parameter  transfer  function  models  as  test  cases.  Specifically,  we  have  been  working  with 
transfer  functions  for  an  Euler- Bernoulli  beam  with  two  kinds  of  damping,  an  Euler-Bernoulli 
beam  with  a  tip  mass,  and  a  Timoshenko  beam,  in  addition  to  a  delay  with  one  stable  pole 
added. 

In  order  to  exactly  solve  the  H°°  optimal  control  problem,  one  must  find  the  inner-outer 
factorization  of  an  H°°  function.  For  the  finite  dimensional  case,  if  one  knows  the  poles  and 
zeros  of  the  transfer  function,  the  factorization  may  be  done  by  inspection;  otherwise  state 
space  algorithms  are  available.  The  latter  approach  is  not  feasible  for  irrational  transfer 
functions,  and  only  a  very  few  simple  examples  may  be  factored  by  inspection.  Therefore, 
we  have  developed  an  algorithm  to  perform  this  factorization  numerically. 

To  justify  this  algorithm,  we  have  proven  two  kinds  of  convergence  results.  The  first 
gives  sufficient  conditions  for  the  inner  factor  of  an  approximant  to  converge  to  the  inner 
factor  of  the  true  function,  whether  the  factorization  is  done  numerically  or  exactly.  This 
result  is  important  since  it  describes  one  set  of  conditions  for  how  well  a  model  can  describe 
the  actual  system  in  order  for  the  i/*”  solution  to  make  sense.  The  second  describes  high 
frequency  behavior  of  the  transfer  function  which  is  sufficient  to  truncate  the  infinite  range 
of  integration  in  the  numerical  procedure.  This  result  is  also  vital  in  describing  the  assump¬ 
tions  which  are  adequate  when  one’s  knowledge  of  the  system  transfer  function  consists  of 
measured  data  and  thus  necessarily  only  contains  information  about  a  finite  frequency  range. 
Finally,  we  have  proven  that  the  H°°  performance  criterion  is  continuous  in  the  inner  factor: 
that  is,  if  a  sequence  of  plant  models  satisfies  the  assumed  convergence  conditions,  the  H°° 
performance  criterion  will  converge  to  that  for  the  true  system.  In  addition  to  validating 
our  computational  techniques,  these  results  answer  the  question  of  well-posedness  of  the  H°° 
control  problem,  providing  conditions  under  which  the  H°°  solution  is  continuous  for  plant 
approximations. 

For  unstable  plants,  a  coprime  factorization  must  be  found  in  order  to  compute  the 
optimal  performance.  In  the  rational  case,  state  space  representations  are  used  with  an 
observer- based  compensator  to  formulate  a  coprime  factorization.  However,  for  irrational 
transfer  functions,  state  space  computation  is  not  an  option.  Instead,  feedback  stabilization 
may  be  used  to  determine  a  coprime  factorization  for  the  plant  from  a  coprime  factorization 
for  a  particular  stabilizing  compensator.  This  method  may  be  applied  experimentally  when 
a  transfer  function  for  the  plant  is  not  explicitly  known  to  arrive  at  a  numerical  coprime 
factorization;  it  may  also  be  used  to  obtain  a  symbolic  coprime  factorization  when  the 
transfer  function  model  is  known.  Analysis  of  the  effect  of  numerical  error  in  this  procedure 
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hu  sdso  been  performed. 

Depending  "upon,  the  problem  formulation,  the  optimal  compensator  may  be  improper, 
which  poses  problems  for  implementation.  We  have  described  a  method  for  approximating 
the  improper  optimal  compensator  with  a  proper  one  and  have  proven  that  one  can  choose 
an  approximation  with  performance  arbitrarily  close  to  optimad.  These  approximation  tech¬ 
niques  have  been  demonstrated  on  a  model  with  an  irrational  outer  factor,  which  could  not 
be  treated  by  previous  methods. 

In  addition,  the  numerical  inner-outer  factorization  is  a  pointwise  computation,  and 
pointwise  knowledge  of  the  inner  factor  leads  to  pointwise  knowledge  of  the  optimal  com¬ 
pensator.  We  would  like  to  have  a  symbolic  representation  of  the  compensator  for  imple¬ 
mentation  purposes;  hence  we  must  interpolate  the  compensator.  Here  we  use  “interpolate” 
to  mean  computing  an  approximation  from  pointwise  information  about  the  true  function. 
We  have  computed  bounds  on  the  permissible  error  in  such  an  approximated  compensator 
to  describe  its  stability  and  performance  properties.  These  bounds  provide  justification  that 
if  we  can  approximate  with  enough  accuracy  (by  increeising  the  number  of  points  at  which 
we  compute  the  optimal  compensator),  we  can  guarantee  the  performance  and  stability  of 
the  resulting  compensator. 
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Appendix  IB 


Research  Summary:  Hong  Yang 

1  Results 

We  summarize  in  the  following  the  main  results  obtained  while  the  author 
was  supported  by  the  grants  mentioned  in  the  title. 

1.  We  observe  that  for  a  general  class  of  SISO  infinite  dimensional  lin¬ 
ear  time  invariant  systems,  including  unstable  systems,  the  resulting 
Hankel  plus  Toeplitz  type  operator  is  a  finite  rank  perturbation  of  the 
Hankel  +  Toeplitz  type  operator  studied  by  Zames  and  Mitter.  We 
show  that  this  finite  riink  perturbation  operator  can  be  computed  ex¬ 
plicitly.  We  also  fill  in  some  missing  points  that  are  required  for  actual 
computation  in  the  Zames-Mitter  work. 

2.  For  a  class  of  MIMO  infinite  dimensional  linear  time  invariant  systems, 
we  give  new  characterizations  for  the  eigenspaces  of  the  corresponding 
Hankel  -f  Toeplitz  type  operator.  We  use  these  new  characterizations  to 
prove  that  the  eigenspaces  are  finite  dimensional  for  this  class  of  Hankel 


+  Toeplitz  type  operator.  Moreover,  we  show  how  to  explicitly  compute 
the  eigenvalues  and  eigenvectors.  We  believe  that  this  is  the  first  time 
that  a  class  of  true  multiple  input-single  output  infinite  dimensional 
linear  time  invariant  systems  is  considered  for  the  H°°  optimal  mixed 
sensitivity  design.  We  also  show  how  to  reduce  the  original  Hankel  -)- 
Toeplitz  type  operator  so  that  its  essential  spectrum  can  be  computed 
explicitly. 

3.  For  computing  the  H”  optimal  mixed  sensitivity  in  the  nonunique 
case,  we  give  new  characterizations  for  certain  subspaces  of  some  Krein 
spaces.  The  new  characterizations  are  used  to  compute  the  H'“  optimal 
mixed  sensitivity  in  some  cases  for  which  it  could  not  be  computed 
before. 

4.  Using  the  characterizations  mentioned  in  3  above,  we  give  some  numer¬ 
ical  schemes  for  computing  the  approximations  of  H*”  optimal  perfor¬ 
mance  and  H"  optimal  pure  and  mixed  sensitivity  for  a  general  class 
of  SISO  systems.  Although  we  have  not  considered  the  convergence  of 
the  algorithms,  the  algorithms  have  been  tested  on  some  infinite  di¬ 
mensional  systems  and  proved  to  be  very  efficient.  A  major  advantage 
of  the  algorithms  is  that  the  complexity  does  not  depend  on  the  di¬ 
mensionalities  of  the  systems  and  weights.  The  algorithms  are  easy  to 
implement. 


5.  We  show  how  to  compute  the  inner-outer  factorization  for  a  class  of 
irrational  H”  matrices.  These  factorizations  are  needed  for  the  H** 
optimal  mixed  sensitivity  design. 

The  fesults  hence  represent  significant  progress  not  only  towards  the  theo¬ 
retical  understanding  of  infinite  dimensional  H”  optimization  problems  but 
also  towards  the  numerical  solution  of  complex  engineering  problems.  We 
note  that  our  approach  is  highly  technical  and  mathematical.  However,  it  is 
also  practical,  as  we  show  with  numerical  examples. 


